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Abstract. We present the results of an ongoing project to use the X-ray observations 
of Type la Supernova Remnants to constrain the physical processes involved in Type la 
Supernova explosions. We use the Tycho Supernova Remnant (SN 1572) as a benchmark 
case, comparing its observed spectrum with models for the X-ray emission from the shocked 
ejecta generated from different kinds of Type la explosions. Both the integrated spectrum of 
Tycho and the spatial distribution of the Fe and Si emission in the remnant are well repro- 
duced by delayed detonation models with stratified ejecta. All the other Type la explosion 
models fail, including well-mixed deflagrations calculated in three dimensions. 
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1. Introduction 

Despite the considerable efforts of the 
last decades, the physical mechanism 
responsible for the explosion of Type 
la Super novae (SNe) still rem ai ns ob - 
scure (Hillebrandt & Niemever, '2000"). 
Traditionally, theoretical calculations of Type 
la SN explosions have been constrained 
mainly through comparison with optical SN 
spectra, but the complexity of the radiative 
transfer calculations involved has only led 



to a limited reduction of the multitude of 
possible explosion mo dels (see, for i n stance 
iBaron et all l2003l) . In iBadenes et all ( l2003h 
and Badene s et alJ (henceforth Paper 

I and Paper II), we introduced a new way to 
explore the physics of Type la SN explosions. 
We have proposed to compare the X-ray 
emission from the shocked ejecta in Type la 
Supernova Remnants (SNRs) with synthetic 
X-ray spectra generated from hydrodynamic 
calculations of the interaction of a grid of Type 
la SN explosion models with the surrounding 
ambient medium (AM), coupled to nonequi- 
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librium ionization (NEI) calculations and an 
X-ray spectral code. 

This comparison is possible because, for 
the first time in history. X-ray observations of 
SNRs with high spatial and spectral resolution 
have become available thanks to Chandra and 
XMM-Newton. Among the SNRs that can be 
observed with sufficient detail, many are still 
young, and their X-ray spectra are dominated 
by the emission from the shocked SN ejecta. 
Several of them are thought to be of Type la 
origin, like the Tycho SNR (SN 1572), SN 
1006 and the Kepler SNR (SN 1604) in our 
Galaxy, or 0509-67.5, DEM L71 and N103B 
in the Magellanic Clouds. Here, we present the 
results from comparison between our models 
and the X-ray spectrum the Tycho SNR. 
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Fig. 1. Spatially integrated XMM-Newton EPIC- 
MOS spectra extracted from two regions (A and B) 
in the western side of the Tycho SNR. Region A en- 
compasses approximately 40% of the SNR, and re- 
gion B is a subset of region A. The most important 
emission lines have been labeled for clarity. 



2. Modeling the X-ray Spectrum of the 
Tycho SNR 

Tycho is one of the brightest X-ray SNRs in 
the sky, and has been observed in great detail 
at all wavelengths from the radio to y-rays. In 
the X-rays, it has been a primary target for all 
the major satellite missions, including ASCA, 
Chandra, and XMM-Newton. Among the con- 
clusions drawn from these observations, we 
highlight the following: 

- The outer edge of emission, which coin- 
cides in the X-rays and in radio observa- 
tions, can be identified with the forward 
shock (FS). 

- The integrated X-ray spectrum of the SNR 
is dominated by strong emission lines from 
Fe, Si, S, Ar, Ca and other elements. 

- The X-ray emission from the shocked AM 
just behind the FS is a featureless con- 
tinuum (Iflwang et al., 2002; Warren et al., 
12005"). This implies that all the X-ray line 
emission has to come from the shocked 
ejecta. 

- The peak of the Fe Ka line emission is at 
lower radii than that of the Fe L or Si Hec 
emission (Hwan g & GotthelA Il997l) . This 
Fe Ka emission comes from material with 
a higher temperature and a lower ioniza- 
tion timescale than the rest of the ejecta 
tHwang et al... 1998.) . 



The XMM-Newton EPIC MOS spectrum 
of the Tycho SNR is presented in Figure [2 
As shown in IWarren et alJ ( l2005l) . the con- 
tinuum emission from the AM is produced 
by a FS whose dynamics is strongly mod- 
ified by cosmic ray (C R) acceleration (see 
iDecourchelle et alll2000l) . This brings the con- 
tact discontinuity between shocked ejecta and 
shocked AM very close to the FS. The re- 
verse shock, on the other hand, is deep into the 
ejecta, and the material close to it is very hot 
(as shown by the morphology of Fe Ka emis- 
sion), so CR acceleration at the reverse shock 
seems unlikely. 

Using these observational results as a start- 
ing point, we attempt to model the X-ray emis- 
sion from the shocked ejecta with the mod- 
els described in Papers I and II. Since all the 
X-ray line emission comes from the shocked 
ejecta, we have used the line flux ratios and line 
centroids determined from the XMM-Newton 
spectrum to reduce the dimensionality of the 
problem and identify the most promising mod- 
els. This is illustrated in Figure |2l where we 
compare the line emission from the Tycho 
SNR with the synthetic spectra for the shocked 
ejecta at a SNR age of f = 430 yr of two 
Type la SN models, a one-dimensional delayed 
detonation (DDTc) and a deflagration calcu- 
lated in three dimensions (B30U). Each panel 
of Fig. 121 corresponds to one of the eight hne 
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Fig. 2. Line flux ratios and line centroids predicted for the delayed detonation model DDTc (left) and the 
three-dimensional deflagration model B30U (right). Values are given as a function of pam (vertical axes), 
and for three different values of fj: p„,i„ (diamonds, solid line), 0.01 (triangles, dotted line), and 0. 1 (squares, 
dashed line). If a particular flux ratio or centroid is missing, this indicates that the relevant line in the 
synthetic spectrum is either very weak or altogether absent. The values determined from the XMM-Newton 
spectrum of the Tycho SNR are plotted as vertical lines, and tolerance ranges are given, with allowed regions 
shaded in grey. 



flux ratios (top two rows of panels) and five 
line centroids (bottom row of panels) that we 
have considered, for a total of thirteen observ- 
able parameters. The observed values are rep- 
resented by vertical lines, with appropriate tol- 
erance ranges and allowed regions shaded in 
grey. For each Type la explosion model, we 
vary the values of pam (the density of the AM 
that is interacting with the ejecta) and (3 (the 
amount of collisionless heating of electrons at 
the reverse shock, defined as the ratio of elec- 
tron to ion postshock specific internal energy 
£e,s/si^s) within ranges that are reasonable for 
the Tycho SNR. See Papers I and II for de- 
tails on the Type la models, definitions of the 
parameters and a complete explanation of the 
simulation scheme. 

Models DDTc and B30U are provided here 
just as representative examples - a more de- 
tailed exploration of the parameter space for 
Type la SN explosions will be given in the 
forthcoming refereed paper As can be seen in 



Fig. 12 model DDTc is much more success- 
ful than B30U. At p^M = 2 ■ lO^-'^g ■ cm^^ 
and 0.01 < /3 < 0.1, model DDTc is capa- 
ble of reproducing twelve of the thirteen ob- 
servable parameters, the only exception be- 
ing the centroid of the Ca Hea line. Model 
B30U, on the other hand, cannot reproduce 
even the fundamental properties of the Fe and 
Si line emission for any combination of pam 
and y6, with important lines like Ar Hea and 
Ca Hea being altogether absent from the syn- 
thetic spectra. When this kind of comparison is 
performed using other models, the conclusion 
is that only one-dimensional delayed detona- 
tions have any hope of reproducing the fun- 
damental characteristics of the line emission 
from the Tycho SNR. All other models fail, 
including one-dimensional deflagrations, sub- 
Chandrasekhar explosions and pulsating de- 
layed detonations. The three-dimensional de- 
flagration models with well-mixed ejecta also 
fail, because the absence of stratification in the 
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ejecta cannot account for the different spec- 
tral properties of Fe and Si (see discussion in 
Section 3 of Paper II). 

Once the most promising values of pam and 
P have been identified through the properties of 
the line emission, the entire synthetic spectra 
can be compared with the observations to eval- 
uate the performance of the models. In Figure 
|3] we plot the synthetic spectra of model DDTc 
(with Pam ^2 - IQ-^'^g- cva^ and (3 = 0.03) 
and model B30U (with pam = 2 ■ 10"-'* g ■ cm"^ 
and p - 0.1), alongside the integrated XMM- 
Newton spectrum. The fit obtained with model 
DDTc is quite impressive, with absolute errors 
comparable to the quality of the atomic data in 
the X-ray spectral codes. It is worth noting that 
the normalization of the flux in the synthetic 
spectra to the observed flux yields an estimate 
for the distance D to the SNR. This 'normal- 
ization distance' D„o„„ is important, because 
the eff'ect of CR acceleration at the FS makes 
it very difficult to constrain D by methods that 
rely on the dynamics of the FS. In the case of 
the DDTc model presented here, Dnom - 2.7 
kpc, a value that is within the 1 .5 < D < 3 kpc 
range proposed bv lSmith et alJ lll99ll) . 

Any successful model for the X-ray emis- 
sion from the Tycho SNR must be able to 
reproduce the spatial distribution of the line 
emission as well as the integrated spectrum. 
In particular, it has to provide an explanation 
for the fact that the Fe Ka line peaks interior 
to both Fe L and Si Hea. In Section 2.4 of 
Paper II we proposed that coUisionless elec- 
tron heating at the reverse shock could pro- 
vide such an explanation, by increasing the 
electron temperature in the innermost shocked 
ejecta and exciting Fe Kff emission close to 
the reverse shock. We illustrate this in Figure 
|4] where we plot the surface brightness pro- 
file of the Si Hea and Fe Ka lines, along with 
that of the 4 to 6 keV continuum emission for 
model DDTc with pam = 2 ■ 10"^"* g ■ cm"^ and 
/3 = 0.03. Qualitatively, these surface bright- 
ness profiles a re similar to those obtaine d by 
XMM-Newton (iDecourchelle et alJ 1200 ll) and 
Chandra (Warren et al., 2005), but a detailed 
quantitative comparison with the data is re- 
quired in order to draw conclusions. This issue 
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Fig. 4. Surface brightness profiles of the Si Hea 
line, 4 to 6 keV continuum and the Fe Kce line for 
model DDTc with p^^ = 2 • 10"-'* g • cm""* and 
P = 0.03. The region of the model represented here 
corresponds mostly to the shocked ejecta: the re- 
verse shock is at a radius of ~ 7.2 ■ 10'* cm and 
the contact discontinuity at ~ 8.2 • 10'** cm. 



will also be addressed in the forthcoming pa- 
per. 

We stress that, even though our analysis 
of the ejecta emission is based on adiabatic 
one-dimensional models, the approximations 
we hav e made are well justified. S orokina et al!| 
1I2OOI have claimed that thermal conduction 
and radiative losses in the SN ejecta need to be 
taken into account for Type la SNRs, but there 
is no observational evidence to support these 
claims. No trace of radiativel y cooled ejecta 
has been found Tycho (see Sm ith et all \\99 1 ). 
and the morphology of the Fe Ka emission can 
only be explained if there is a temperature gra- 
dient in the shocked ejecta, which is incompat- 
ible with efficient thermal conduction. 



3. Conclusions 

In this contribution, we have summarized the 
main results of our ongoing efforts to use 
the X-ray emission from the Tycho SNR to 
constrain the physics of Type la SN explo- 
sions. These results are very promising: we 
have found that only one class of models, one- 
dimensional delayed detonations, is capable 
of reproducing the fundamental properties of 
the ejecta emission. Among the delayed det- 
onations, the preferred model is DDTc, with 
Ek - 1.16-10^' erg and 0.8 solar masses 
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Fig. 3. Comparison of the synthetic spectra with the best combination of p/^^ (in g • cm"^) and fi, for models 
DDTc (left) and B30U (right) with the XMM-Newton observation of t he Tycho SNR. T he dotted line in the 
DDTc plot represents the power law continuum with index 2.72 (after lFink et al^ll994^ added to model the 
contribution from the shocked AM emission. The continuum from the shocked ejecta in model B30U is so 
high that the fit did not admit an additional power law component. 



of Ni. This model provides a good approx- 
imation to the integrated X-ray spectrum for 
PAM = 2 ■ 10-24 g . cm-'' and /? = 0.03. Other 
constraints, like the surface brightness profile 
of the most relevant lines, the normalization 
distance D„orm, and the apparent radii of the 
contact discontinuity and the reverse shock, are 
also in agreement with the observations (see 
[Warren et al., 2005). More details will be pro- 
vided in the forthcoming refereed paper. 

It is remarkable that the sophisticated 
three-dimensional deflagration models with 
well -mixed eject a published by several groups 
(see lBravo et al.L 12005. for a review) fail to re- 
produce the fundamental properties of the X- 
ray emission from Tycho, in particular the fact 
that Fe and Si have different spectral prop- 
erties. This indicates that some degree of el- 
emental stratification must exist in Type la 
SN ejecta. The ultimate physical cause of this 
stratification remains an open issue, and even 
though some explanations have been proposed, 
no self-consistent model for Type la explosions 
exists yet. 
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